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Abstract 

A framework for developing new approximate electronic structure methods is presented, in which 
the correlation energy of a many-electron system in the ground state is computed as in the single- 
reference second-order many-body perturbation theory but with the reference one-body Hamilto- 
nian modified by additional (generally non-local) potentials which are universal functionals of the 
Hartree-Fock density. The existence of such functionals that reproduce the exact correlation en- 
ergy justifies the search for approximate models which may overcome some important deficiencies 
of the traditional density-functional methods and also perform better than the usual low-order 
wavefunction methods. 
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I. INTRODUCTION 



The Hohenberg-Kohn theorems l| have provided a formal justification for the develop- 
ment of a number of approximate density functionals for electronic structure calculations 
within the Kohn-Sham scheme Q|. Such density-functional methods are now widely used 
in theoretical chemistry for computing ground state properties of molecules due to their 
low cost and scalability. Nevertheless, there is still a number of deficiencies that the tra- 
ditional approximate functionals cannot overcome. Self-interaction error Q is one such 
example, another important problem is the inability to describe the dispersion (van der 
Waals) interaction |4|. On the other hand, the simplest correlated wavefunction method - 
the M0ller-Plesset second-order (MP2) many-body perturbation theory j^J - does not suffer 
from these deficiencies, its cost, though higher, is still acceptable for many interesting ap- 
plications. There are cases, however, where this simple approximation shows its weakness - 
for strongly-correlated systems the correlation energy can be significantly overestimated, an- 
other difficulty appears when the underlying Hartree-Fock reference experiences an artificial 
spin- or simmetry-breaking. 

In this work we propose a framework for developing new approximated methods of elec- 
tronic structure calculations which combine the strengths of both the low-order many-body 
perturbation theory and the local density functional approximations. 

II. THEORY 

Starting from the fact that the exact total energy of a many-electron system in the ground 
state is a functional of the Hartree-Fock density [6j , we propose here a special form of the 
correlation energy functional for which simple approximate models may be developed. 

The single-determinat wavefunction $ of the Hartree-Fock theory assumed normalized 

(*{*) = 1 (1) 
makes the expectation value of the energy 

E = ($\H\$) (2) 

stationary, which is equivalent to the conditions 

($t\Hm=0, (3) 



where singly-substituted determinants $" are introduced. The indices k, . . . (a,b,c, . . .) 
label the occupied (virtual) one-electron wavefunctions, whereas k, A, fi, v will refer to any 
of them. 

In terms of one- and two-electron integrals 

H; = f 0,(r) (-|V 2 + „(r)) 0,(r) d 3 r, (4) 

p^ _ pAf p^A 

^KfJ, -'"'K/i - rl 'K/i' \°) 

R% = J 0.(r)^(r')ir - r'rVA(r)i(r') d 3 r d 3 r', (6) 

where the wavefunctions K (r) are assumed to be real and their spin dependence is not show 
explicitly for simplicity, the energy ([2]) and the conditions take the form 

Eo = HI + iJ3g, (7) 

J? = 0, (8) 

with the Fock matrix defined by 

f; = h; + b%. (9) 

The summation over all repeated indices is assumed throughout. 

We now introduce a first-order correlated wavefunction ^ as a linear combination of all 
doubly substituted determinants 

*! = ±7#<&#, (10) 
from which the (exact) correlation energy will be computed as 

£ c =(*il#l$> (11) 

or in terms of two-electron integrals 

E c = \t$B$. (12) 
The coefficients are determined by linear equations 

fa^cb i fb^ac fk^ab rk^ab i r>a6 n / i o\ 

Jc + Jc^j ~ Ji T kj - J j T ik + R ij = (13) 

which are formally derived from many-body perturbation theory with the one-electron ma- 
trix elements of the zeroth-order Hamiltonian in the form 



./y n • «■!. (14) 

fa = F'a + <, (15) 
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u\ = J ^(r) Mo (r;[p ])^(r)d 3 r, (16) 
u\ = J a (r) Mv (r;[ Po ])0 6 (r)d 3 r. (17) 

The two different "correction potentials" u Q and u v to be added to the one-electron Fock 
operator are universal functionals of the Hartree-Fock density 

p (r) = &(r)&(r) (18) 



and have the property that the resulting correlation energy (1111) is equal to the exact one. 
It is obvious that the correction potentials thus defined are not unique. The simplest form 
would be a constant shift (independent of r) for either all occupied or all virtual energy 
levels: u Q (r) = u[po], u v = or vice versa, but in that case the density functional u[po] 
should be highly non-local and very hard to model without violating the size-consistency 
of an approximate method. One would be interested in the models where these potentials 
were local in nature. 

For the full-CI wavefunction of the form 

ci = $ + + \T* ^ + h T $®tk + ■■■ ( 19 ) 

the correlation energy can be expressed as ([12]) in terms of the doubles coefficients Tf b 
alone, so it would be natural to require the equations ( IT3|) to generate the values rf b which 
approximate T® b in ( TO?]) . Having only n (n + l)/2 + n v (n v + l)/2 — 1 matrix elements u\ 
and u b a as parameters (n Q and n v are the dimensions of the occupied and virtual subspaces, 
the latter can be in principle infinite) it is not possible, in general, to reproduce all n (n — 
l)n v (n v — l)/4 nontrivial values of T?-. The differences between the reduced one-electron 
quantities 

Q J i = T$r? k - T%T% (20) 
Q^ = r^-T^ (21) 

can be brought to zero by an appropriate choice of ul and u b a in ( fl3l) . in that case, however, 
the correlation energy ( Till may differ from the exact one. The minimization of either 

q = QlQl + Q b a Q b a (22) 

or 

q = QiHlQk + QaRltQt (23) 
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with respect to u\ and u h a with the constraint of reproducing the full-CI correlation energy 
by ( I12p leads to the equations which uniquely define the correction potentials in either case. 
A more simple function 

q = (t$ - T$) {(F?5 bd + F b d 8 ac )8 ik 8 fl - SMtfSn + F 1 ^)) (iff - (24) 

can also be used as a measure of the deviation between the doubles coefficients. 

The existence of such (not necessarily multiplicative) potentials u a (r; [p ]) and w v (r; [po]) 
as functionals of the Hartree-Fock density follows from the fact that there is a one-to- 
one mapping between the external potential v(r) and the Hartree-Fock density Po( r )> the 
relationship can be schematically described as p =>■ v =>• ^ci u ,u v . 

In analogy to the traditional density-functional theory the development of approximate 
functionals of this new type can start from the consideration of the uniform electron gas as 
the simplest model system. A local density approximation of the form 

Uo (r) = (po(r)) a (p (r)) , (25) 
u v (r) = U ^(p (r))(l + a(po(r))), (26) 

can be proposed, where the function u neg (p) is chosen as to reproduce, within this method, 
the exact correlation energy of the uniform electron gas, and the function a(p) can be 
parametrized in some way, the simplest choices being the constants a = or a = -1. 

III. CONCLUSIONS 

The new form of the correlation energy functional of the Hartree-Fock density proposed 
here is based on a physically meaningful description of the electron correlation using the 
two-particle equations (|T3|) . Approximate methods of this type should be able to describe 
the dispersion interaction in a natural way and should not contain self-interaction error. It 
remains to be seen if the simplest local density approximation ( 125]) and (1261) will show higher 
accuracy in molecular applications than the standard MP2 method. Such comparison will 
become possible as soon as the function u neg (p) will be available from numerical calcula- 
tions of the uniform electron gas, which we are currently investigating. More sophisticated 
density-functional models of our corrections potentials u Q (r; [p ]) and u v (r; [p ]) can also be 
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developed, an empirical parameterization [7| being the simplest choice. 
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